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ABSTRACT 



Asteroid families, traditionally defined as clusters of objects in orbital parameter space, often have 

l> ■ 

£T) • distinctive optical colors. We show that the separation of family members from background inter- 
im 

lopers can be improved with the aid of SDSS colors as a qualifier for family membership. Based on 

00 " 

an ~88,000 object subset of the Sloan Digital Sky Survey Moving Object Catalog 4 with available 

>: 

proper orbital elements, we define 37 statistically robust asteroid families with at least 100 members 
(12 families have over 1000 members) using a simple Gaussian distribution model in both orbital 
and color space. The interloper rejection rate based on colors is typically ~10% for a given orbital 
family definition, with four families that can be reliably isolated only with the aid of colors. About 
50% of all objects in this data set belong to families, and this fraction varies from about 35% for 
objects brighter than an H magnitude of 13 and rises to 60% for objects fainter than this. The 
fraction of C-type objects in families decreases with increasing H magnitude for H > 13, while the 
fraction of S-type objects above this limit remains effectively constant. This suggests that S-type 
objects require a shorter timescale for equilibrating the background and family size distributions 
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via collisional processing. The size distribution varies significantly among families, and is typically 
different from size distributions for background populations. The size distributions for 15 families 
display a well-defined change of slope and can be modeled as a "broken" double power-law. Such 
"broken" size distributions are twice as likely for S-type familes than for C-type families (73% vs. 
36%), and are dominated by dynamically old families. The remaining families with size distributions 
that can be modeled as a single power law are dominated by young families (<1 Gyr). When size 
distribution requires a double power-law model, the two slopes are correlated and are steeper for 
S-type families. No such slope-color correlation is discernible for families whose size distribution fol- 
lows a single power law. For several very populous families, we find that the size distribution varies 
with the distance from the core in orbital-color space, such that small objects are more prevalent 
in the family outskirts. This "size sorting" is consistent with predictions based on the Yarkovsky 
effect. 
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1 Introduction 

The size distribution of asteroids is one of most significant observational constraints on their history 
and is considered to be the "planetary holy grail" (Jedicke & Metcalfe 1998, and references therein). 
It is also one of the hardest quantities to determine observationally because of strong selection effects. 
Recently, Ivezic et al. (2001, hereafter 101) determined the asteroid size distribution to a sub-km 
limit using multi-color photometry obtained by the Sloan Digital Sky Survey (York et al. 2000; 
hereafter SDSS). Here we extend their work by using an updated (4 th ) version of the SDSS Moving 
Object Catalog (Ivezic et al. 2002a, hereafter I02a). 

The main goal of this paper is to study size distributions of asteroid families. Asteroid dynamical 
families are groups of asteroids in orbital element space (Gradie, Chapman & Williams 1979, Gradie, 
Chapman & Tedesco 1989, Valsecchi et al. 1989). This clustering was first discovered by Hirayama 
(1918, for a review see Binzel 1994), who also proposed that families may be the remnants of 
parent bodies that broke into fragments. About half of all known asteroids are believed to belong to 
families; Zappala et al. 1995 (hereafter Z95), applying a hierarchical clustering method to a sample 
of 12,487 asteroids, finds over 30 families. Using the same method and a larger sample of ~106,000 
objects, Nesvorny et al. (2005, hereafter N05) identify ~50 statistically robust asteroid families. 

The size distributions of asteroid families encode information about their formation and evolution, 
and constrain the properties of the families' parent bodies (e.g., Marzari, Farinella & Davis 1999; 
Tanga et al. 1999; Campo Bagatin & Petit 2001; Michel et al. 2002; de Elia & Brunini 2007; Durda 
et al. 2007; and references therein). Motivated by this rich information content, as well as the 
availability of new massive datasets, here we address the following questions 

(1) What is the fraction of objects associated with families? 

(2) Do objects that are not associated with families show any heliocentric color gradient? 

(3) Do objects that are not associated with families have uniform size distribution independent of 
heliocentric distance? 



(4) Do objects associated with families have a different size distribution than those that are not 
in families? 

(5) Do different families have similar size distributions? 

(6) Is the size distribution related to family color and age? 



These questions have already been addressed numerous times (e.g. Mikami & Ishida 1990; Cellino, 
Zappala & Farinella 1991; Marzari, Davis & Vanzani 1995; Z95; Morbidelli et al. 2003; N05). The 
main advantages of the size distribution analysis presented here, when compared to previous work, 
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The large sample size: we use a set of ~88,000 objects for which both SDSS colors and proper 
orbital elements computed by Milani & Knezevic (1994) are available 

Simple and well-understood selection effects: the SDSS sample is >90% complete without a strong 
dependence on magnitude (Juric et al. 2002, hereafter J02) 

Improved faint limit: the sample of known objects listed in the latest ASTORB le from Jan- 
uary 2008 to which SDSS observations are matched is now essentially complete to r ~ 19.5 
(corresponding to H ~ 17 in the inner belt, and to H ~ 15 in the outer belt) 
Improved family definitions due to color constraints (rejection of interlopers and separation of 
families overlapping in orbital space) 
Improved accuracy of absolute magnitudes derived using SDSS photometry, as described below. 



The SDSS asteroid data are described in Section 2, and in Section 3 we describe a novel method for 
defining asteroid families using both orbital parameters and colors. Analysis of the size distribution 
for families and background objects is presented in Section 4, and we summarize our results in 
Section 5. 



2 SDSS Observations of Moving Objects 

2. 1 An Overview of SDSS 

The SDSS is a digital photometric and spectroscopic survey which will cover about one quarter of 
the Celestial Sphere in the North Galactic cap, and produce a smaller area (~300 deg 2 ) but much 
deeper survey in the Southern Galactic hemisphere (Stoughton et al. 2002; Abazajian et al. 2003, 
2004, 2005; Adelman-McCarthy et al. 2006). SDSS is using a dedicated 2.5m telescope (Gunn et 
al. 2006) to provide homogeneous and deep (r < 22.5) photometry in five bandpasses (Fukugita 
et al. 1996; Gunn et al. 1998; Smith et al. 2002; Hogg et al. 2001; Tucker et al. 2006) repeatable 
to 0.02 mag (root-mean-square scatter, hereafter rms, for sources not limited by photon statistics, 
Ivezic et al. 2003) and with a zeropoint uncertainty of ~0. 02-0. 03 mag (Ivezic et al. 2004). The 
flux densities of detected objects are measured almost simultaneously in five bands (u, g, r, i, and 
z) with effective wavelengths of 3540 A, 4760 A, 6280 A, 7690 A, and 9250 A. The large survey 
sky coverage will result in photometric measurements for well over 100 million stars and a similar 
number of galaxied_|. The completeness of SDSS catalogs for point sources is ~99.3% at the bright 
end and drops to 95% at magnitudes of 22.1, 22.4, 22.1, 21.2, and 20.3 in u, g, r, % and z, respectively. 
Astrometric positions are accurate to better than 0.1 arcsec per coordinate (rms) for sources with 
r < 20.5 (Pier et al. 2003), and the morphological information from the images allows reliable star- 
galaxy separation tor ~ 21.5 (Lupton et al. 2002, Scranton et al. 2002). A compendium of other 



technical details about SDSS can be found on the SDSS web site (http://www.sdss.org), which also 
provides interface for the public data access. 



1 The recent Data Release 6 lists photometric data for 287 million unique objects observed in 9583 deg 2 



of sky; see http://www.sdss.org/dr6/. 



2.2 SDSS Moving Object Catalog 

The SDSS, although primarily designed for observations of extragalactic objects, is significantly 
contributing to studies of the solar system objects because asteroids in the imaging survey must 
be explicitly detected and measured to avoid contamination of the samples of extragalactic objects 
selected for spectroscopy. Preliminary analysis of SDSS commissioning data by 101 showed that 
SDSS will increase the number of asteroids with accurate five-color photometry by more than 
two orders of magnitude, and to a limit about five magnitudes fainter (seven magnitudes when the 
completeness limits are compared) than previous multi-color surveys (e.g. The Eight Color Asteroid 
Survey, Zellner, Tholen & Tedesco 1985). For example, a comparison of SDSS sample with the Small 
Main-Belt Asteroid Spectroscopic Survey (Xu et al. 1995; Bus & Binzel 2002ab) is discussed in detail 
by N05. 

SDSS Moving Object Catalog__| (hereafter SDSS MOC) is a public, value-added catalog of SDSS 
asteroid observations (I02a). It includes all unresolved objects brighter than r = 21.5 and with 
observed angular velocity in the 0.05-0.5 deg/day interval. In addition to providing SDSS astro- 
metric and photometric measurements, all observations are matched to known objects listed in the 
ASTORB file (Bowell 2001), and to a database of proper orbital elements (Milani 1999; Milani 
& Knezevic 1994)), as described in detail by J02. J02 determined that the catalog completeness 
(number of moving objects detected by the software that are included in the catalog, divided by the 
total number of moving objects recorded in the images) is about 95%, and its contamination rate 
is about 6% (the number of entries that are not moving objects, but rather instrumental artifacts). 
The most recent SDSS MOC 4th data release contains measurements for 471,000 moving objects. A 
subset of 220,000 observations were matched to 104,000 unique objects listed in the ASTORB file 
(Bowell 2001). The large sample size increase between the first and fourth release of SDSS MOC 
is summarized in Figure 1. The object counts in both releases are well described by the following 



http://www.sdss.org/dr6/products/value_added/index.html 



function (101) 
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where x = r — re, a — (ki + A; 2 )/2, b — (ki — k 2 )/2, with k\ and k<i the asymptotic slopes of log(n) 
vs. r relations. This function smoothly changes its slope around r^rc, and we find best-fit values 
re = 18.5, k\ = 0.6 and ki = 0.2. The normalization constant, n , is 7.1 times larger for SDSS 
MOC 4 than for the first release. In addition to this sample size increase, the faint completeness 
limit for objects listed in ASTORB also improved by a about a magnitude, to r ~ 19.5 (the number 
of unique ASTORB objects increased from ~11,000 to ~100,000). Above this completeness limit, 
the SDSS MOC lists color information for ~33% of objects listed in ASTORB. 

The quality of SDSS MOC data was discussed in detail by 101 and J02, including a determination of 
the size and color distributions for main-belt asteroids. An analysis of the strong correlation between 
colors and the main-belt asteroid dynamical families was presented by Ivezic et al. (2002b, hereafter 
I02b). Jedicke et al. (2004) reported a correlation between the family dynamical age and its mean 
color for S-type families, and proposed that it is due to space weathering effects. This correlation was 
further discussed and extended to C-type families by N05. Multiple SDSS observations of objects 
with known orbital parameters can be accurately linked, and thus SDSS MOC also contains rich 
information about asteroid color variability, discussed in detail by Szabo et al. (2004) and Szabo 
and Kiss (2008). 



2.3 Errors in H magnitudes listed in the ASTORB file 

As pointed out by J02, there is a large systematic discrepancy between the absolute magnitudes 
listed in ASTORB file and values implied by SDSS measurements. The latter are computed as 

H CO rr = H ASTORB + V — cV, (2) 



where Hastorb is the ASTORB value, cV is the apparent magnitude in Johnson system com- 
puted from information listed in ASTORB as described in J02, and V is the observed magnitude 
synthesized from SDSS g and r magnitudes (SDSS MOC entries 47, 42, and 32, respectively). 

This discrepancy persists in the 4 th release of SDSS MOC, as illustrated in Figure 2. The mean 
difference between H measured by SDSS and the values from ASTORB is 0.23 mag, and the root- 
mean-scatter is 0.30 mag. The best-fit shown in Figure 2 implies that uncertainty of H corr is about 
0.16 mag, with a negligible systematic error (the latter is expected to be about 0.02-0.03 mag due 
to uncertainties in absolute photometric calibration of SDSS imaging data; see section 2.1). It is 
likely that this uncertainty is dominated by magnitude variation due to rotation. The magnitude 
offset of 0.33 mag for ~70% of measurements implied by the best fit could be due to measurements 
reported by LINEAR. A similar magnitude offset at the faint end is a known problem in LINEAR 
calibration, and is currently being addressed with the aid of new calibration catalogs (J.S. Stuart, 
priv. comm.). 

Since the random error in H is twice as large as for H corr , we adopt H corr in the remainder of this 
work. For a detailed analysis of this magnitude offset problem! 3 1, we refer the reader to J02. 



3 The Asteroid Families in SDSS MOC 

The contrast between dynamical asteroid families and the background population is especially strong 
in the space defined by proper orbital elements. These elements are nearly invariants of motion and 
are thus well suiteouLlfor discovering objects with common dynamical history (Valsecchi et al. 1989, 
Milani & Knezevic 1992). 



IAU Commission 15 has formed a Task Group on Asteroid Magnitudes to address this problem, see 



http://www.casleo.gov.ar/cl5-wg/index-tgh.html 
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The current asteroid motion is usually described by osculating orbital elements which vary with time 
due to perturbations caused by planets, and are thus less suitable for studying dynamical families. 
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The value of SDSS photometric data becomes particularly evident when exploring the correlation 
between colors and orbital parameters for main-belt asteroids. I02b demonstrated that asteroid 
dynamical families, defined as clusters in orbital element space, also strongly segregate in color 
space. We use the technique developed by I02b to visualize this correlation for ~45,000 unique main- 
belt asteroids with H corr < 16 listed in SDSS MOC 4 (Figures 3-4). The asteroid color distribution 
in SDSS bands shown in Figure 3, and its comparison to traditional taxonomic classifications, is 
quantitatively discussed by 101 and N05. 

A striking feature of Figure 4 is the color homogeneity and distinctiveness displayed by asteroid 
families. In particular, the three major asteroid families (Eos, Koronis, and Themis), together with 
the Vesta family, correspond to taxonomic classes K, S, C, and V, respectively (following Burbine 
et al. 2001, we assume that the Eos family is associated with the K class). Their distinctive op- 
tical colors indicate that the variations in surface chemical composition within a family are much 
smaller than the compositional differences between families, and vividly demonstrate that asteroids 
belonging to a particular family have a common origin. 

3.1 A Method for Defining Families Using Orbits and Colors 

Traditionally, the asteroid families are defined as clusters of objects in orbital element space. The 
most popular methods for cluster definition are the hierarchical clustering and the wavelet analysis 
(Z95, N05). Given the strong color segregation of families, it is plausible that SDSS colors can be 
used to improve the orbital family definitions and minimize the mixing of candidate family members 
and background population. 

The SDSS colors used to construct Figures 3-4 are the i — z color and the so-called a* color, defined 
in 101 as 

a* = 0.89(g-r) + 0A5(r-i)-0.57 (3) 

The a* color is the first principal component of the asteroid color distribution in the SDSS r — 



% vs. g — r color-color diagram (for transformations between the SDSS and Johnson system see 
Ivezic et al. 2007). Similar principal component analysis was also performed Roig and Gil-Hutton 
(2006), who considered the distribution of taxonomic classes (especially V-type asteroids) in SDSS 
principle components by comparing directly to spectroscopic data, and by N05, whose two principal 
components are well correlated with the a* and i — z colors (we find that a* = 0.49 PC\ — 0.16 
reproduces the measured a* values with an rms of 0.026 mag for objects with r < 18). The principal 
colors derived by N05 include the u band, which becomes noisy at the faint end. Given that the 
completeness of the known object catalog (ASTORB) reaches a faint limit where this noise becomes 
important, we use the a* and i — z colors to parametrize the asteroid color distribution. Therefore, 
the family search is performed in a five- dimensional space defined by these two colors and the proper 
semi-major axis, sine of the inclination angle and eccentricity. 

There are numerous techniques that could be used to search for clustering in a multi-dimensional 
space (e.g. Z95; N05; Carruba & Michtchenko 2007). They differ in the level of supervision and 
assumptions about underlying data distribution. Critical assumptions are the distribution shape for 
each coordinate, their correlations, and the number of independent components. We utilize three 
different methods, one supervised and two fully automatic. The automatic unsupervised methods 
are based on the publicly available code FASTMLXI 5 1 by A. Moore and a custom- written code based 
on Bayesian non-parameteric techniques (Ferguson 1973; Antoniak 1974). 

In the supervised method (1) families are manually identified and modeled as orthogonal (i.e. aligned 
with the coordinate axes) Gaussian distributions in orbital and color space. The two unsupervised 
methods (2 and 3) also assume Gaussian distributions, but the orientation of individual Gaussians 
is arbitrary, and the optimal number of families is determined by the code itself. All three methods 
produce fairly similar results and here we describe only the supervised method (1), and use its 
results in subsequent analysis. The two unsupervised methods produce generally similar results for 
the objects associated with families, but tend to overclassify the background into numerous (50-60) 
small families, and their details and results are not presented thoroughly in this paper. 



See http://www.cs.cmu.edu/~psand 
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We select from the SDSS MOC 4 the first observation of all objects identified in ASTORB, and 
for which proper orbital elements are also available, resulting in 87,610 objects. Among these, 
there are 45,502 objects with H corr < 16. We split the main sample into three subsets using semi- 
major axis ranges defined by the major Kirkwood gaps (see Figure 4): inner (a < 2.50), middle 
(2.50 < a < 2.82) and outer (a > 2.82) belt. For each subset, we produce the e vs. sin{i) diagrams 
color-coded analogously to Figure 4, and use them to obtain preliminary identification of asteroid 
families in both orbital and color space. Approximate rectangular bounds are assigned to these 
visually identified families, from which median (centroid) and standard deviation, a, for the three 
orbital elements are estimated. Using these estimates, for each asteroid we compute distance in 
orbital space from a given family centroid as 



Dorm = yjd* + d 2 e + d* (4) 



\flcentroid ^object) 



u a 


o~ a 


d e = 


\&centroid ^object) 


cr e 


di = 


[^centroid ^object) 


°i 



where 

(n . .,-n ,. .\ 

(5) 
(6) 
(7) 

Histograms in D or .b it were used to determine a preliminary value of A<>6, the maximum orbital 
distance from a family centroid for an object to be ascribed family membership. These initial Ao6 
are determined from the differential D or ut distribution as the position on the first local minimum 
(the further rise of counts with increasing D or b it is due to the background objects and other fam- 
ilies). The object distribution in the D or bit vs. a* and D or ut vs. {% — z) diagrams was used to first 
define approximate rectangular bounds for each family, and then to compute the color centroid and 
standard deviation in a* and (i — z) for each candidate family. In order to use color as a family 
discriminator, we define analogously to orbital elements 



D co i or = Jdj + d\ (8) 
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where 



(7" /j* 

i _ \ centroid object) /q\ 

j _ l,* ^centroid ''^object) /-. «\ 



We use histograms in D co i or to define Ac i, the maximum color distance from a family centroid for 
an object to be ascribed family membership. Figure 8 illustrates the D or bu and D co i or histograms 
for the Vesta and Baptistina families, and the bottom panels show the distributions of family and 
background objects in the D or bit vs. Dcoior plane. 

In cases of families which formed from the disruption of a differentiated parent body the color 
distribution might not provide a well defined morphology. However, we did not find any case where 
a subset of objects selected using D or bit did not result in one or two well defined color distributions. 
Nevertheless, it is possible that a small fraction of objects could be rejected from a family due to 
different colors than the majority of other members. 

All objects that have both D or ut < &-Orb and D co i or < ^Coi are then considered to be a family 
member. With a given estimate of family populations, this procedure is iterated and all parameters 
are refined. It typically takes one to two iterations to converge. All "converged" families are removed 
from the sample and the process was repeated until there were no family candidates with more than 
100 members. This condition is the result of the requirement that the statistical errors for the slope 
of absolute magnitude distribution of the families are smaller (typically 0.01-0.02) than plausible 
systematic errors (0.03-0.04), as discussed below. 

Using this procedure, we found 37 families which account for 46% of all objects. Their defining 
parameters are listed in Table 1. Additional three candidate families that had fewer than 100 
members in the last iteration were discarded (see last three entries in Table 2). The family names 
were determined by comparison with Z95 and N05, and, when no corresponding family was found, 
by searching for the lowest-numbered asteroid in the Milani & Knezevic (1994) catalog of proper 
orbital elements. In a small number of cases, it is possible that the "name-giving" object has a color 
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that is inconsistent with the majority of objects in the family. We ignored such cases and retained 
nomenclature from Z95 and N05 in order to ease the comparison. That is, the defining properties 
of families are the orbital and color parameters listed in Table 1, rather than their names. 

The separation of the main-belt asteroids into families and the background objects is illustrated in 
Figures 5 and 6. While there is some residual structure in the background, both in color and orbital 
space, it is much less prominent than the structure displayed by identified families. 

3.2 The Comparison of Resulting Families with Previous Work 

Important questions about the quality of family associations derived here, that are relevant for 
subsequent analysis, are 

• The contamination: are all families robust? 

• The completeness: are any families missed? 

• What is the impact of variations in adopted family definitions on the resulting family properties 
such as the number of members and their size distribution? 

In this section we address the first two questions, and discuss the third one below (§4.3). 

It is unlikely that the derived sample of 37 families contains any spurious family because of the 
conservative requirement that a family must include at least 100 members. Indeed, we have rejected 
three familes that are likely real because thay had fewer candidate members (5, 90 and 46). The very 
narrow color distribution of all selected families provides another argument for their robustness. 

Our method of iterative removal of identified families and the simultaneous use of colors and orbital 
parameters to search for remaining families appears very robust when compared to two automated 
methods also employed. Nevertheless, given that all three methods from this work assume gaussian 
distributions of colors and orbital parameters, it is prudent to compare our family list with families 
obtained by other means, such as the hierarchical clustering method. We have cross-correlated the 
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list of 37 families determined here with the list of 41 families obtained by N05 using the hierarchical 
clustering method. N05 based their study on a larger sample of objects with proper orbital elements 
(~106,000 vs. ~88,000 analyzed here; note that the latter sample extends to ~1.5 mag fainter 
flux limit but is smaller because it includes only objects observed by SDSS), and did not place a 
requirement on the minimum number of objects per family. Therefore, it is plausible that families 
missed by our selected method may be present in their list. 

Out of 41 families from the N05 list, 27 are listed in Table 2. This is encouraging level of agreement 
given the significant difference in applied methodology. We examined in detail each of the fourteen 
N05 families missing from our list and searched for them in the sample of background objects. We 
did not find any candidate family that included more than 100 members, though most appear to 
be real clusters. 

Among the ten families from our list that we could not identify in the N05 list, three were detected 
by at least one method discussed by Z95, and thus are likely real (Euterpe, Teutonia, and Henan). 
Of the remaining seven, the recognition of four families was greatly aided by color information 
(Baptistina from Flora, Mitidika and Juno, Lydia and Padua, and McCuskey from Nysa-Polana) . 
It is likely that the remaining three families were not detected by N05 because they have steep 
absolute magnitude distributions and thus only a small number of members were present in the 
(older) version of catalog used by N05. For example, among the 3405 objects in the Teutonia 
family, only 37 have H corr < 14. 

We further compared our list of families to those presented in Mothe-Diniz et al. (2005), who used 
spectroscopic measurements to probe asteroid family structure. Of their 21 nominal families, all 
but three (Renate, Hoffmeister, and Meliboea) can be matched to families detected here. Of the 
numerous smaller (many with fewer than 100 members) "clumps" they identify, 10 correspond 
to families listed here, while eight families (Teutonia, Mitidika, Euterpe, Andree, Lydia, Ursula, 
Lyxaohua and Theobaldia) present in our list do not appear in Mothe-Diniz et al. (2005). They 
resolved the remaining family in our data set, Flora, into a number of smaller clumps that merged 
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into a single family at high cutoff velocities. 

A good example of the separation of dynamically mixed families using SDSS colors is provided 
by the small family Baptistina which is "buried" within the Flora family. Figures 7 and 8 (right 
panels) illustrate how different a* color distributions enable the identification of ~5% of the objects 
nominally assigned Flora family ((a*) = 0.13) membership by their orbital parameters as being 
members of the Baptistina family ((a*) = —0.04). While initially puzzled why a similar color-aided 
search by N05 did not yield any additional families, we have found that all color-separated families 
extracted here are dominated by faint objects and thus may not have been present in sufficiently 
large numbers in the older catalog. We conclude that all the families discussed here are robustly 
detected, and that it is very unlikely that we missed any family with more than 100 hundred 
members. It is, however, possible that the background is composed of numerous families dominated 
by small objects that are not discernible with the presently available catalog. Hence, the fraction of 
~ 50% of objects associated with families is only a lower limit (this fraction is a function of object 
size, as discussed below). 

Figure 11 shows the color dependence of the family and background populations on semi-major axis. 
We note that the median a* color for the background population becomes bluer as the semi-major 
axis increases in the same way as the median color for the family population. 

4 The Size Distribution of Asteroid Populations 

The known object catalog is complete to r ~ 19.5; above this limit an ASTORB entry is found for 
practically every SDSS moving object. Depending on the distance and orientation of the observed 
object, this apparent magnitude limit corresponds to a completenes limit ranging from H ~ 17 in 
the inner belt to H ~ 15 in the outer belt. Brighter than these limits, selection function is essentially 
equal to 1 for the purposes of this work (SDSS managed to observe only about 1/3 of all ASTORB 
objects, but this is essentially a random selection without an impact on derived absolute magnitude 
distribution of individual families). These simple selection effects allow us to derive robust family 
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size distributions to very small size limits. The only other study of family size distributions that 
approached the same size limit is that of Morbidelli et al. (2003), who had to introduce an ad 
hoc indirect correction for selection effects in the known object catalog (which is supported by our 
analysis, as discussed below). 

The transformation between the asteroid absolute magnitude, H, and its effective diameter, D, 
requires the knowledge of the absolute visual albedo pv, 

H = 18.1 - 2.5 log(-) -5log(D/lkm). (11) 

The absolute albedo is not known for the overwhelming majority of objects in our sample. However, 
the albedo is known to be strongly correlated with colors (Zellner 1979; Shoemaker et al. 1979; 
Delbo 2004); for example, the C-like asteroids (a* < 0) have a median albedo of 0.04 and the S-like 
asteroids have a median albedo of 0.14. Given that the color variations within a family are small, it 
seems plausible that the albedo variation within a family is also small (this is supported by the data 
compiled by Tedesco, Cellino & Zappala 2005; see their Table 7). With this assumption, the shapes 
of the absolute magnitude distribution and the distribution of log(-D) are the same. Hereafter, we 
will interchangebly use "the absolute magnitude distribution" and "size distribution", where the 
latter implies the distribution of log(D). For simplicity, in the remainder of analysis we only use 
the differential distributions. 

If the differential absolute magnitude distribution, n(H) = AN/ AH, can be described by 

log(n) = Const. + aH, (12) 

and the albedos of objects within a given family or population are similar, then it follows from eq. 11 
that the differential size distribution can be described as n(D) oc D~ q , with the size distribution 
index 

q = 5a + l. (13) 
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While the absolute magnitude distributions derived here often cannot be described by a single 
power-law, eq. 13 is still useful for locally relating the slope of the H distribution to the slope 
of implied (differential) D distribution. For example, a model based on an equilibrium cascade in 
self-similar collisions developed by Dohnanyi (1969) predicts q = 3.5 and a = 0.5. 



4-1 The Comparison of Size Distributions for Families and Background 

We compare the size distributions for the family population and for the background, separately for 
the three regions defined by semi- major axis. The differential absolute magnitude distributions are 
shown in Figure 9. To aid the comparison of different panels, we plot for reference the differential 
distribution derived from the cumulative distribution reported by 101 

IQax 

n(r) = n — — , (14) 

where x = H corr — H c , a = (ki + fc 2 )/2, b = (&i — k 2 )/2, with H c = 15.5, k\ = 0.65 and k 2 = 0.25 
(Table 4 in 101). 

101 were able to fit this functional form because their sample extended to a ~1.5 mag fainter 
H limit (H corr ~ 17.5) than the sample discussed here. Given this sample difference, for each H 
distribution shown in Figure 9 we instead fit a "broken" power law: a separate power-law fit for the 
bright and faint end. While this procedure is expected to yield a shallower slope at the faint end 
than the above 101 fit, it is preferred here because it "decouples" the bright and faint ends. The 
separation of the bright and faint ends was attempted in H steps of 0.5 mag, and the value that 
minimizes the resulting \ 2 was adopted as the best fit. The statistical errors for the best-fit slopes 
are typically 0.01-0.02, but it is likely that their uncertainty is perhaps a factor of two or so larger 
due to systematic effects (see §4.3 below). In a few cases, the best fit is consistent with a single power 
law. The best-fit power-law parameters for differential absolute magnitude distributions shown in 
Figure 9 are listed in Table 3. 
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The data and best fits shown in the top three panels in Figure 9 demonstrate that the absolute 
magnitude distributions are not identical: the outer main-belt shows a flatter distribution, and the 
inner belt shows a steeper distribution than the middle belt region for objects with H corr < 14. This 
is in conflict with the 101 finding that the size distribution appears universal throughout the belt. 
However, here we analyze a sample about seven times larger; the statistical errors at the bright end 
for the 101 sample were too large to detect this effect (see their Figures 21 and 22). Nevertheless, 
the 101 size distribution remains valid when the whole belt is treated together because the counts 
underprediction of their fit in the outer belt is compensated by its overprediction in the inner belt. 

The separation of populations into families and background (the middle and bottom rows in Fig- 
ure 9) shows that the flattening of H corr distribution as the semi-major axis increases is valid for 
each subpopulation separately. Objects associated with families always show the flattening at the 
faint end, while the background populations admit a single power-law fit in the middle and outer 
belt. 

Due to different H corr distributions for family and background populations, the fraction of objects 
associated with families is a function of H corr . The top left panel in Figure 10 shows this dependence 
separately for blue (dominated by the C taxonomic type) and red (dominated by the S type) 
subsets. For both subsets, the fraction of objects in families significantly increases from ~20% to 
~50% between H corr = 9 and H corr = 11. The two color-selected subsets show different behavior for 
Hcorr > 11: for blue subset the fraction of objects in families decreases from ~50% to ~30%, while it 
stays constant at the ~60% level for red families. Since blue families typically have larger semi-major 
axis than red families, it is possible that this decrease in family membership is due to increasing 
color rejection at the faint end. However, the remaining two panels in Figure 10 demonstrate that 
this is not the case because the color rejection rate is both fairly independent of H corr , and too 
small to account for the observed decrease of blue family membership. 

The dominance of background objects for H corr < 13 is consistent with the background population 
having a significantly shallower size distribution for large objects than the families. The falloff of blue 
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family fraction and the slow climb in color rejection rate toward large values of H corr in Figure 10 
confirms that the size distributions for blue families are shallower for values of H corr > 13. Because 
the red family fraction is effectively flat to our detection limit it appears that the red family and 
background populations have identical size distributions for objects with H corr > 13. Morbidelli et 
al (2003) suggest that the background population is composed of many small families which formed 
from small-diameter objects and as such should have steep size distributions for objects larger than 1 
km, producing a background population with an initial size distribution (for objects >1 km) steeper 
than that for families which formed from the breakup of larger objects. Differences between the size 
distributions of these two populations should eventually disappear through collisional processing. 
Because the blue family and background populations appear to have significantly different size 
distributions while the red populations' size distributions appear identical (for H corr > 13), we infer 
that the time for equilibrating the background and family size distributions is longer for blue objects 
than red. This difference in the equilibration time may be due either to differences in asteroid internal 
structure and material properties between taxonomic classes or due to environmental variations in 
collisional processing rates, as red objects are more prevalent in the inner belt and blue objects in 
the outer belt. 

We note that the background population in the outer belt shows a curious excess of large objects 
(H corr < 11.5) compared to best power-law fit (Figure 9, bottom right). We have inspected the 
orbital parameter and color distributions for 58 objects with 10 < H corr < 11 and found that they 
are not associated with any identified family, nor generally clustered. 



4-2 The Comparison of Size Distributions for Individual Families 

The inspection of differential H corr distributions for the 37 families identified here shows that many, 
but not all, display a clear change of slope such as seen for family populations in Figure 9. We have 
attempted a "broken" power-law fit for all families. When the two best-fit slopes differ by less than 
0.05, we enforce a single power-law fit. This procedure yields 22 families described by a single power 
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law and 15 families with a robust detection of the slope change. Their best-fit parameters are listed 
in Tables 4 and 5, respectively, and a few examples of measured H distributions and best fits are 
shown in Figure 12. 

For families whose absolute magnitude distributions are described by a single power law, the median 
best-fit power-law slope is 0.56, with a standard deviation of 0.16 (determined form inter-quartile 
range). This scatter is significantly larger than the measurement errors and indicate that families 
do not have a universal size distribution. Similarly, for families with a best-fit "broken" power law, 
the medians and standard deviations for the "bright" and "faint" slopes are (0.66, 0.24) and (0.32, 
0.15), respectively (again note the significant scatter relative to the measurement errors), with the 
median H corr where the slope changes of 14.2 (D ~6 km for py = 0.1). We discuss correlations of 
these best-fit parameters with the family color and age in §4.4. 

4-3 Systematic Deviations in Size Distribution due to Variations in Family Definitions 

Before proceeding with the analysis of correlations between size distributions and other family 
properties such as color and age, we analyze the systematic deviations in size distribution due to 
variations in family definitions. For example, the color constraints may result in a size-dependent 
incompleteness because of the increased photometric noise at the faint end. Similarly, the assumption 
of gaussian distributions for orbital parameters and colors may result in incomplete families due 
to extended halos, as pointed out by N05. This effect may also induce size- dependant systematics 
because small objects are scattered over a larger region of orbital space, as shown below. 

The Vesta family offers a good test case because of its unique color distribution (which is due to 
the influence of 1 /im absorption feature on the measured % — z color). The top panel in Figure 13 
compares the H corr distributions for the adopted Vesta family and for a less constraining orbital 
cut defined simply by 0.06 < sin(i) < 0.16 and e < 0.16, that yields 30% more candidate members. 
Apart from this overall shift in the normalization, the resulting distributions have statistically indis- 
tinguishable shapes. The middle panel compares the adopted family and a much less constraining 

20 



color cut, a* > (i.e. no constraint on the i — z color), that yields 50% more objects. Again, the 
slope of the two distributions are indistinguishable. 

We detect a significant difference, however, when we split the adopted family in the "core" and 
"outskirt" parts using D orbit < 1 and 1.75 < D or t> it < 2.75 (see eq. 4). As the bottom panel 
in Figure 13 shows, the "outskirt" subsample has a steeper H corr distribution than the "core" 
subsample. The best-fit power-law slopes in the 14 < H corr < 16 region are 0.45 and 0.59 for the 
"core" and "outskirt" subsample, respectively. We note that despite this slope difference, the change 
of slope between the bright (0.89) and faint end is robustly detected. 

Another method to see the same "size sorting" effect is to inspect the dependence of D orbit on H corr . 
We find that the median D or ut for objects in the Vesta family increases from 1.0 to 1.5 as H corr 
increases from 14 to 17. 

This "size sorting" effect is not a peculiar property of the Vesta family as it is seen for a large fraction 
of families. It is caused by an increased scatter in all three orbital parameters as H corr increases. 
This is not surprising as the velocity field of the fragments produced in the disruption of an asteroid 
family's parent body may have been size- dependent. For most families the sorting is dominated by 
the increased dispersion in the semi-major axis. One of the most striking examples, the Eos family, 
is shown in Figure 14. As discussed by N05, this increase of dispersion as size decreases can be also 
be explained as the drift induced by the Yarkovsky effect (see also Vokrouhlicky 1999; and Bottke 
et al. 2001). 

4-4 Correlations between Size Distributions and Family Color and Age 

We analyze the correlations between the best-fit size distribution parameters listed in Tables 4 and 
5, and family color and age. The age, when available, is taken from the compilation by N05. 

The dependence of the power-law index on the mean a* color for families described by a single 
power law is shown in the top left panel in Figure 15. The mean and standard deviation for 14 blue 
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families are (0.55, 0.13), and for 8 red families are (0.65, 0.19). These differences are not statistically 
significant. Within each color-selected subsample (blue vs. red, i.e. a* < vs a* > 0), there is no 
discernible correlation between the slope and color. 

Families that require a "broken" power law fit are twice as likely for red families (a* > 0) dominated 
by S type asteroids) than for blue families dominated by C type asteroids (73% vs. 36%). As 
illustrated in Figure 15, the size distributions are systematically steeper for S type families, and 
the "bright" and "faint" end slopes appear to be correlated. The median values of the "bright" and 
"faint" end slopes are (0.57, 0.18) for blue families, and (0.79, 0.39) for red families. 

For a subset of families that have available age estimates, we find that families with "broken" power 
law size distributions are dominated by old families, while those that admit a single power law are 
dominated by young families, with the age separation boundary at ~1 Gyr. We note that the size 
distribution was used for some of age estimates compiled by N05, so this conclusion may be a bit of 
circular reasoning, though the majority of age estimates are derived independently of the observed 
size distribution. 

The correlations between the mean color and family age reported by Jedicke et al. (2004) and 
N05 are reproduced when using the a* color for families discussed here. Figure 16 illustrates a 
good agreement with the analytic fits to the observed correlations obtained by N05, and further 
demonstrates the correlation between the observed size distributions and age. 



5 Discussion and Conclusions 

We have used a large sample of asteroids (~88,000) for which both orbital elements and SDSS colors 
are available to derive improved membership for 37 asteroid families. The addition of colors typically 
rejects about 10% of all dynamically identified candidate members due to mismatched colors. Four 
families can be reliably isolated only with the aid of colors. About 50% of objects in this data set 
belong to families, with this fraction representing a lower limit due to a conservative requirement 
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that a candidate family must include at least 100 members. The resulting family definitions are in 
good agreement with previous work (e.g. Z95, N05) and all the discrepancies are well understood. 
Although SDSS has observed only about 1/3 of all known asteroids, it is remarkable that the 
sample discussed here provides color information for more than an order of magnitude more objects 
associated with families than analyzed in the published literature. 

This data set enables the determination of absolute magnitude (size) distributions for individual 
families to a very faint limit without a need to account for complex selection effects. We verify that 
size distribution varies significantly among families, and is typically different from size distributions 
for background populations. Consequently, the asteroid size distribution cannot be described by a 
universal function that is valid throughout the main belt (e.g. Jedicke & Metcalfe 1998, Ivezic et al. 
2001, and reference therein). This finding will have an influence on conclusions derived from mod- 
eling the size distribution under this assumption (e.g. Bottke et al. 2005, and references therein). 
In particular, it is not clear how to interpret a detailed dependence of the critical specific energy 
(energy per unit mass required to fragment an asteroid and disperse the fragments to infinity) on as- 
teroid size derived from such models, when the starting observational constraint on size distribution 
is an average over multiple families with significantly varying size distributions. 

We show that for objects with H corr < 13, the background population dominates (family fraction 
decreases toward lower H corr , indicating a shallower size distribution for large objects), while for 
objects with H corr > 13, the red family fraction remains effectively constant to our completeness 
limit while the blue family fraction falls off. This indicates that the time to collisionally equilibrate 
the family and background populations (see e.g., Morbidelli et al. 2003) is shorter for red objects 
than blue. 

The size distributions for 15 families display a well-defined change of slope and can be modeled as 
a "broken" double power-law. The first evidence for this effect and a discussion of its significance 
are presented by Morbidelli et al. (2003). Using a data set with much simpler correction for the 
observational selection effects, we confirm their result in a statistically more robust way. We also 
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find such "broken" size distributions are twice as likely for S-type familes than for C-type families 
(73% vs. 36%), and are dominated by dynamically old families. The remaining families with size 
distributions that can be modeled as a single power law are dominated by young families (<1 Gyr). 

The eight largest families all show a change of size distribution slope to much smaller values at 
the faint end (see Table 5). This result has a direct consequence when prediciting the number of 
very small objects (D ~ 1 km). In particular, it could explain why the Statistical Asteroid Model 
developed by Tedesco, Cellino & Zappala (2005) predicts too many objects: the data presented here 
are inconsistent with the SAM assumptions for the number of objects in its most populous families 
such as Eunomia and Themis. 

We find that when size distribution requires a double power-law model, the two slopes are correlated 
and are steeper for S-type families. No such slope-color correlation is discernible for families whose 
size distribution follows a single power law. While beyond the scope of this work, the modeling of 
such correlations may shed light on the internal structure and material properties of asteroids. 

For several very populous families, we find that the size distribution varies with the distance from 
the core in orbital-color space, such that small objects are more prevalent in the family outskirts. 
As discussed by N05 (and references therein), this "size sorting" is consistent with predictions based 
on the Yarkovsky/YORP effect. 

While these results provide significant new observational constraints for the properties of main-belt 
asteroids, very soon the observations will further improve. The upcoming large-scale sky surveys, 
such as Pan-STARRS (Kaiser et al. 2002) and LSST (Tyson 2002), will obtain even more impressive 
samples, both in size, diversity of measurements and their accuracy. For example, LSST will scan 
the whole observable sky every three nights in two bands to a 5cr depth equivalent to V = 24.7. 
These data will enable much improved analysis due to several factors 

• Due to hundreds of observations, the orbits will be determined directly, instead of relying on 
external data, resulting in a sample about 30-40 times larger than discussed here 

• The effective faint limit will be extended by about 5 magnitudes, correspoding to ten times smaller 
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size limit (diameters of several hundred meters) 

• Due to many photometric observations obtained with the same well-calibrated system, the un- 
certainties in absolute magnitudes will be an order of magnitude smaller 

• The addition of the y band (at ~1 /im) will improve the color classification due to better sensitivity 
to the ~1 jum absorption feature present in spectra of many asteroids. 

These new data will undoubtely reinvigorate both observational and theoretical studies of main-belt 
asteroids. 
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Fig. 1. An illustration of the improvements in sample size between the first and the fourth release of SDSS 
Moving Object Catalog. Symbols with (statistical) error bars show differential counts for moving objects 
listed in the first release from 2002 (dots: all objects detected by SDSS; circles: identified in ASTORB file) 
and the fourth release from 2008 (triangles: all SDSS; squares: in ASTORB). Note that, in addition to 
a sample size increase of about a factor of 7, the faint completeness limit for objects listed in ASTORB 
also improved by a about a magnitude (the number of unique ASTORB objects increased from ~11,000 
to ~100,000). The dashed lines show a double-power law fit described in text, with the dlog(N)/dr slope 
changing from 0.60 at the bright end to 0.20 at the faint end. For illustration, the dot-dashed lines shows 
a single power-law with a slope of 0.60. 
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Fig. 2. A comparison of asteroid absolute magnitude, H, inferred from SDSS measurements, and the value 
listed in ASTORB file for ~133,000 observations of about 64,000 unique objects observed at phase angles be- 
tween 3 and 15 degrees. The histogram shows the data distribution (AH = Hsdss—Hastorb = H corr —H, 
see eq. 2) and the dot-dashed line is a best fit. The best fit is a linear combination of three gaussians: two 
(which simulate asteroid variability) are centered on 0.02, have widths of 0.08 and 0.20 mag, and have 
relative normalizations of 13% and 18%, respectively. Their sum is shown by the dashed line centered on 
AH = 0.02. The third gaussian (which accounts for a large number of objects with bad photometry) has 
a width of 0.28 mag and is shown by the dashed line centered on AH = 0.33. That is, about 69% of H 
measurements listed in ASTORB file are systematically too bright by 0.33 mag. 
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Fig. 3. A plot of the color distribution in a* and i — z of 45,087 unique objects listed in both the SDSS 
MOC 4 and ASTORB file, and that have H corr < 16. The approximate boundaries of three spectral classes 
are marked, and used in labeling family type. The color-coding scheme defined here is used in figures 4-6. 




Fig. 4. A plot of the proper a vs. sin(i) for the same objects as shown in Figure 3. The color of each dot is 
representative of the object's color measured by SDSS, according to the color scheme defined in Figure 3. 
The three main regions of the belt, defined by strong Kirkwood gaps, are marked. 
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Fig. 5. Illustration of the decomposition of the main-belt asteroid population into families and background 
objects in proper a vs. sin(i) (left panels) and proper a vs. e (right panels). The top panels show all 
(background and family) objects in the data subset. The two middle panels show objects from 37 identified 
families, and the bottom two panels show the background population. 



35 




Fig. 6. Analogous to Figures 5, except that the top three panels show the e vs sin{i) distribution for the 
three main regions defined by strong Kirkwood gaps (a < 2.5 left, 2.5 < a < 2.82 middle, 2.82 < a < 3.5 
right; see Figure 4). The middle row shows family members (with several families of note labeled), and the 
bottom row shows the background population. For a high-resolution version of this figure with complete 



labeling, see http://www.astro.washington.edu/ivezic/sdssmoc/sdssmoc.html 
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Fig. 7. An illustration of color differences for families with practically identical orbital parameter distri- 
butions. The dashed histogram shows the a* color distribution for 6,164 candidate members of the Flora 
family. The solid histogram shows the a* color distribution for 310 candidate members of the Baptistina 
family, which is easily separated from the Flora family thanks to the SDSS color information. 
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Fig. 8. Top Left: D^^i histogram for the Vesta family and surrounding objects (defined from the Vesta 
family centroid). Top Right: Equivalent to top left plot, but for the dynamically buried Baptistina family. 
The dotted line represents the D or bi t histogram without any color constraints, which climbs smoothly to 
very high numbers of objects because of the inclusion of Flora family objects. The solid line represents 
the D or i,i t histogram with the color constraints applied, which displays a stronger clustering signature. 
Vertical dashed line represents Ao r b cutoff values selected for these families. Middle Left: D co i or histogram 
for objects that met the D or f,i t criteria for the Vesta family. Middle Right: D co i or histogram for objects in a 
preliminary orbital definition of the Baptistina family, showing strong color distinction from the background 
Flora objects. Vertical dashed line represents Acoi cutoff values selected for these families. Bottom Left 
and Right: D or ut vs. D co i or for Vesta and Baptistina families, respectively. Dashed box defines Aob and 

Acoi boundaries (objects inside are assigned family membership). 
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Fig. 9. The differential absolute magnitude distributions corresponding to panels in Figure 6 are shown as 
symbols with (Poisson) error bars. The solid line shows arbitrarily renormalized best-fit distribution from 
101. The two dashed lines show the best-fit "broken" power law: a separate power-law fit for the bright and 
faint end. In some cases, the two lines are indistinguishable. The best-fit parameters are listed in Table 3. 
The two arrows show the best-fit break magnitude (left) and the adopted completeness limit (right). 
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Fig. 10. The family-to-background population ratios and effects of color selection on family populations. 
The top left panel shows the fraction of objects in families to the total population as a function of H corr 
magnitude, with the solid histogram representing red objects (a* > 0) and the dotted histogram repre- 
senting blue objects (a* < 0). The top right panel compares the fraction of objects in families to the total 
population as a function of magnitude when colors are used as a constraint on family membership (solid 
histogram) and when they are not (dotted histogram). The bottom panel shows the rejected fraction due 
to color constraints as a function of magnitude. 
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Fig. 11. Median a* color as a function of semi-major axis for several populations. Solid histogram represents 
the family population, dashed line represents the background population, and the dotted line represents 
the median of the family a* color centroids (unweighted by family size) for each Main Belt region (inner, 
mid, and outer). 
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Fig. 12. Analogous to Figure 9, except that the absolute magnitude distributions for selected asteroid 
families are shown. The first three panels (from top left to bottom right) show examples of families that 
follow a single power-law magnitude distribution, and the remaining six panels show magnitude distribu- 
tions for families that require a double power-law fit. The best-fit parameters are listed in Tables 4 and 5, 
respectively. 
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Adopted Vesta vs. orb. relaxed Vesta 




Fig. 13. The dependence of absolute magnitude distribution on family definition for Vesta family. The 
top panel compares the magnitude distribution obtained for adopted family definition (open circles) to 
that obtained for a much less constraining cut on orbital elements (closed squares). The middle panel 
compares the adopted distribution (open circles) to that obtained for a much less constraining color cut 
(closed squares). In both cases the shapes of the magnitude distributions are similar, with only significant 
change in the number of selected candidate members. The bottom panel separates the adopted population 
into the orbital "core" region (open circles) and "outskirt" region (closed squares). Note that the latter 
distribution is steeper (i.e. the small members are more prevalent in the "outskirt" region). 
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effect. Note that this family is intersected by several mean motion resonances with Jupiter (e.g. 7:3 at 
a p ~ 2.95). 
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Fig. 15. A summary of relationships between the family median color and parameters of the best-fit 
power-law magnitude distributions. The top left panel shows the slope of the best-fit power-law as a 
function of color a* (see text for definition) for 25 families that follow a single power-law magnitude 
distribution. The top right panel shows the "bright" slope of the best-fit double power-law as a function of 
color a* for 12 families that follow a double power-law magnitude distribution, and the bottom left panel 
is analogous plot for the "faint" slope. Note that blue (C type) families have much shallower magnitude 
distributions than redder families. The bottom right panel demonstrates the strong correlation between 
the "bright" and "faint" slopes, that seems independent of color. The errors of alpha are dominated by 
systematics. As per discussion in Section 4.3, we estimate these systematic errors to be about 0.05-0.1. 
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Fig. 16. The correlation between the mean a* color for families defined here and their age taken from 
N05. Familes whose absolute magnitude distribution can be described by a single power law are shown 
by circles, and those that require a broken power law as squares. Note that the former are dominated by 
young families (<1 Gyr), and the latter by old families. The two lines are the best-fit color-age relation 
from N05, converted using a* = 0.49 PC\ — 0.16, where PC\ is the first principal color component derived 
by N05 (note that N05 used several very young familes not discussed here to constrain the slopes of plotted 
relations). 
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3.00 


2.50 


2.91 


3.60 


Themis 1 


3.126 


0.071 


0.024 


0.008 


0.153 


0.019 


-0.120 


0.037 


0.010 


0.059 


2.50 


2.50 


2.91 


3.60 


Hygiea 1 


3.144 


0.066 


0.090 


0.008 


0.137 


0.019 


-0.110 


0.037 


0.010 


0.082 


2.50 


2.50 


2.91 


3.60 


Lixiaohua 1 


3.149 


0.015 


0.177 


0.004 


0.197 


0.005 


-0.070 


0.037 


0.080 


0.067 


3.00 


2.50 


2.91 


3.60 


Ursula 


3.156 


0.060 


0.279 


0.012 


0.090 


0.028 


-0.070 


0.044 


0.060 


0.067 


2.50 


3.00 


2.91 


3.60 


Veritas 1 


3.168 


0.007 


0.160 


0.003 


0.063 


0.004 


-0.080 


0.030 


0.050 


0.052 


3.00 


3.50 


2.91 


3.60 


Thcobalda 


3.170 


0.014 


0.247 


0.006 


0.251 


0.011 


-0.160 


0.037 


0.040 


0.089 


2.75 


3.50 


3.00 


3.50 



Notes: 

a the family name (the lowest-numbered member), sorted by semi-major axis 

b the mean proper semi-major axis (AU) and its gaussian dispersion (a) adopted for this family 

c the mean sine of proper inclination and its gaussian dispersion adopted for this family 

d the mean proper orbital eccentricity and its gaussian dispersion adopted for this family 

e the mean a* color (see text for definition) and its gaussian dispersion adopted for this family 

' the mean i — z color and its gaussian dispersion adopted for this family 

9 the maximum deviation in orbital space from the adopted mean values, in units of adopted dispersions (see text) 

h the maximum deviation in color space from the adopted mean values, in units of adopted dispersions (see text) 

1 the minimum and maximum semi-major axis adopted for this family 

1 family matched to N05 



Table 1 

Adopted family definitions in the orbital and color space. 
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Family 


N b 


a c p 


;d 


c p 


a/ 


i - z° 


Flora 


6164 


2.28 


0.08 


0.13 


0.13 


-0.05 


Baptistina 


310 


2.28 


0.10 


0.14 


-0.04 


0.03 


Vesta 


3793 


2.35 


0.11 


0.10 


0.12 


-0.32 


McCuskey 


1043 


2.36 


0.05 


0.15 


-0.12 


0.01 


Erigone 


307 


2.37 


0.09 


0.21 


-0.09 


0.05 


Euterpe 


387 


2.38 


0.02 


0.18 


0.09 


-0.03 


Nysa-Polana 


2928 


2.39 


0.04 


0.18 


0.13 


-0.04 


Andree 


649 


2.40 


0.08 


0.17 


0.13 


-0.02 


Massalia 


730 


2.41 


0.03 


0.16 


0.07 


-0.04 


Teutonia 


3405 


2.57 


0.07 


0.16 


0.12 


-0.04 


Rafita 


225 


2.59 


0.13 


0.17 


0.08 


-0.04 


Maria 


1315 


2.61 


0.25 


0.09 


0.12 


-0.02 


Mitidika 


698 


2.61 


0.22 


0.24 


-0.11 


0.03 


Eunomia 


2995 


2.63 


0.23 


0.15 


0.13 


-0.03 


Misa 


185 


2.65 


0.04 


0.18 


-0.08 


0.05 


Adeona 


428 


2.66 


0.20 


0.17 


-0.11 


0.05 


Juno 


354 


2.66 


0.23 


0.23 


0.07 


-0.03 


Aeolea 


172 


2.66 


0.07 


0.19 


-0.09 


0.04 


Henan 


624 


2.67 


0.04 


0.06 


0.11 


-0.02 


Nemesis 


129 


2.73 


0.09 


0.09 


-0.09 


0.02 


Lydia 


598 


2.74 


0.11 


0.03 


0.16 


0.01 


Padua 


442 


2.75 


0.09 


0.04 


-0.05 


0.05 


Merxia 


252 


2.75 


0.09 


0.14 


0.08 


-0.07 


Gefion 


914 


2.76 


0.16 


0.13 


0.10 


-0.03 


Chloris 


121 


2.76 


0.15 


0.25 


-0.06 


0.05 


Agnia 


1106 


2.76 


0.06 


0.07 


0.08 


-0.01 


Dora 


248 


2.79 


0.14 


0.20 


-0.12 


0.04 


Brasilia 


127 


2.85 


0.26 


0.12 


-0.05 


0.03 


Koronis 


1267 


2.90 


0.04 


0.05 


0.09 


-0.02 


Eos 


4367 


3.04 


0.18 


0.07 


0.04 


0.03 


Tirela 


411 


3.12 


0.29 


0.20 


0.15 


0.08 


Themis 


1073 


3.13 


0.02 


0.15 


-0.11 


0.01 


Hygica 


1076 


3.15 


0.09 


0.13 


-0.11 


0.01 


Lixiaolma 


150 


3.15 


0.18 


0.20 


-0.07 


0.05 


Ursula 


644 


3.15 


0.28 


0.09 


-0.06 


0.06 


Veritas 


250 


3.17 


0.16 


0.06 


-0.08 


0.05 


Theobalda 


100 


3.17 


0.25 


0.25 


-0.15 


0.01 


NA h 


90 


2.87 


0.11 


0.05 


0.17 


0.04 


NA h 


5 


2.77 


0.08 


0.09 


-0.08 


0.14 


NA h 


46 


2.78 


0.06 


0.07 


0.00 


-0.25 



Notes: 

a the family name (the lowest-numbered member) 

° the number of objects in SDSS MOC 4 associated with this family 

c the median proper semi-major axis (AU) 

d the median sin of proper inclination 

e the median proper orbital eccentricity 

f the median a* color (see text for definition) 

9 the median i — z color 
h rejected clumps with less than 100 members 

Table 2 

The median orbital parameters and SDSS colors for detected families and rejected clumps (last three). 
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a range" 


N 6 


H c - 


■LJC 

rL max 


H| 


«! 


4 


2.00 - 2.50 


30,702 


11.0 


16.5 


14.0 


0.76 


0.46 


family 


16,309 


11.0 


16.0 


13.5 


0.79 


0.53 


background 


14,393 


11.0 


16.0 


13.5 


0.69 


0.46 


2.50 - 2.82 


32,500 


11.0 


16.0 


13.5 


0.73 


0.42 


family 


14,261 


11.0 


16.0 


13.5 


0.67 


0.44 


background 


18,239 


11.0 


15.5 


13.5 


0.54 


0.57 


2.82 - 3.60 


24,367 


12.0 


15.5 


13.5 


0.56 


0.40 


family 


9,547 


12.0 


15.0 


13.5 


0.57 


0.37 


background 


14,820 


12.0 


15.0 


13.5 


0.56 


0.52 



Notes: 

a the range of proper semi-major axis for defining the inner, middle and outer main belt (AU); 

b the number of objects in each subsample; the first line corresponds to the full sample, and the following two to subsamples classified 

as families and background, respectively. The total number of objects is 87,569. 

c the minimum and maximum H magnitude used in fitting the H distribution 

d the best-fit "break" H magnitude (see text) 

e the "bright" H distribution slope 

f the "faint" H distribution slope 

Table 3 

Best-fit parameters for counts shown in Fig. 9 



Family" 


N b 


H c b 


Hf 


"1 


age-' 


Baptistina 


310 


14.0 


16.0 


0.49 


— 


McCuskey 


1043 


12.5 


16.0 


0.49 


— 


Erigone 


307 


14.0 


16.0 


0.59 


— 


Euterpe 


387 


13.5 


16.0 


0.52 


— 


Andree 


649 


14.0 


16.0 


0.70 


— 


Massalia 


730 


14.5 


16.0 


0.97 


0.3 ± 0.1 


Rafita 


225 


14.0 


16.0 


0.35 


1.5 ± 0.5 


Mitidika 


698 


12.5 


15.5 


0.61 


— 


Misa 


185 


13.0 


15.5 


0.47 


0.5 ± 0.2 


Juno 


354 


14.5 


16.0 


0.63 


— 


Aeolea 


172 


14.0 


15.5 


0.85 


— 


Nemesis 


129 


14.0 


16.0 


0.70 


0.2 ± 0.1 


Lydia 


598 


14.0 


16.0 


0.81 


— 


Padua 


442 


14.0 


16.0 


0.50 


— 


Merxia 


252 


14.5 


15.5 


0.71 


0.5 ± 0.2 


Chloris 


121 


13.0 


15.5 


0.35 


0.7 ± 0.4 


Agnia 


1106 


13.0 


15.5 


0.52 


0.2 ± 0.1 


Brasilia 


127 


13.5 


15.5 


0.71 


0.05 ± 0.04 


Lixiaohua 


150 


13.0 


15.0 


0.56 


0.3 ± 0.2 


Ursula 


644 


11.0 


15.0 


0.47 


— 


Theobalda 


100 


13.0 


15.5 


0.44 


— 


Veritas 


250 


12.0 


15.0 


0.50 


8.3 ± 0.5 Myr 



Notes: 

a the family name (the lowest-numbered member) 

" the number of objects in SDSS MOC 4 associated with this family 

c bright H magnitude limit used for fitting 

d faint H magnitude limit used for fitting 

e the best-fit power-law index for H distribution 

e the family age in Gyr (except for Veritas in Myr), taken from Ncsvorny et al. (2005), when available 

Table 4 

Best-fit H distribution parameters for families described by a single power law. 
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Family" 


N b 


H t 


Hf 


h b 


4 


«2 


age h 


Flora 


6164 


11.0 


16.0 


14.0 


0.66 


0.40 


1.1 ± 0.5 


Vesta 


3793 


12.5 


16.0 


14.8 


0.89 


0.50 


— 


NysaPolana 


2928 


13.5 


16.5 


15.0 


0.80 


0.39 


— 


Teutonia 


3405 


12.5 


15.5 


14.0 


0.94 


0.59 


— 


Maria 


1315 


11.0 


15.5 


14.5 


0.53 


0.38 


3.0 ± 1.0 


Eunomia 


2995 


11.0 


15.5 


13.5 


0.76 


0.28 


2.5 ± 0.5 


Adeona 


428 


12.0 


15.5 


14.2 


0.57 


0.29 


0.7 ± 0.5 


Henan 


624 


14.0 


16.0 


15.1 


0.79 


0.46 


— 


Gefion 


914 


12.0 


15.5 


14.2 


0.86 


0.26 


1.2 ± 0.4 


Dora 


248 


13.0 


16.0 


14.8 


0.37 


0.18 


0.5 ± 0.2 


Koronis 


1267 


12.0 


15.0 


13.5 


0.55 


0.26 


2.5 ± 1.0 


Eos 


4367 


10.5 


15.0 


13.9 


0.52 


0.32 


2.0 ± 0.5 


Tirela 


411 


13.2 


15.0 


14.0 


1.04 


0.62 


— 


Themis 


1073 


11.0 


15.0 


13.4 


0.46 


0.10 


2.5 ± 1.0 


Hygiea 


1076 


12.7 


15.0 


14.4 


0.62 


0.18 


2.0 ± 1.0 



Notes: 

a the family name (the lowest-numbered member) 

the number of objects in SDSS MOC 4 associated with this family 

c bright H magnitude limit used for fitting 

d faint H magnitude limit used for fitting 

e the "break" H magnitude used for fitting 

' the best-fit power-law index for H distribution at the bright end 

9 the best-fit power-law index for H distribution at the faint end 

h the family age in Gyr, taken from Nesvorny et al. (2005), when available 

Table 5 

Best-fit H distribution parameters for families described by a "broken" power law. 
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